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THE CALCULATION OF ANTENNA RADIATION PATTERNS 
BY A VECTOR THEORY USING DIGITAL COMPUTERS 


Richard F. Schmidt 
Advanced Development Division 


ABSTRACT 


This document Continues with the application of the field equations to the 
solution of antenna problems begun in X-525-68-201. The electric and magnetic 
fields at the reflector boundary after reflection are related to the incident fields 
via the boundary conditions. It is shown that the evaluation of the scattered 
radiation pattern can be accomplished in an efficient manner by converting 
electric field terms to their magnetic equivalents and utilizing only a simple 
magnetic polarization vector. The history of the wave emanating from the 
prime-feed source is next written into the formulation explicitly, but without 
introducing a specific illumination function. A generally-applicable technique 
for accommodating feed translation and rotation with six degrees of freedom is 
then developed for numerical integration of the field equations. An appendix 
containing the subroutines for surface normal and differential area of the set of 
conic surfaces (paraboloid, hyperboloid, ellipsoid, sphere, and cone) is included 
so that this document provides the means for evaluating a large class of antenna 
configurations. 
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GLOSSARY OF NOTATION 


S3nnbol 

Meaning 

V 

vector operator del 

r.R 

radial distances, subscripted 

0(r) 

scalar wave without time dependence 

n 

normal to surface or line, a vector 

S, s 

surface, subscripted 

C 

contour 

OJ 

angular frequency 

x,y,z;x', y', z'; x", y”, z" 

Cartesian coordinates for reference frames 

k 

wave number 

E, iff 

electric and magnetic field vectors, subscripted 

■ 

electrical conductivity, subscripted 

. 

magnetic permeability, subscripted 

e 

inductive capacity, subscripted 

D, B 

dielectric displacement and magnetic flux density 

V 

surface charge density 

K 

sheet current, subscripted 

dT 

differential length, a vector 

j 

a constant equal to in complex variable theory 


unit vectors in the r, ^, 0 directions, spherical 
coordinate system 
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Symbol 

a, /3, y 

Pe 

G(0, <I>) 
e. , h. 

cr, 5 

F 

a, c 
^1 


GLOSSARY OF NOTATION (Continued) 

Meaning 

Euler angles of the rigid body transforniation 
feed translation vector 
feed gain function 

electric and magnetic polarization vectors 

radial and angular variables of the cylindrical 
coordinate system 

focal length, and certain defined functions 

paiameters of hyperboloids, ellipsoids, spheres, 
cones 

curvature of a curve 



THE CALCULATION OF ANTENNA RADIATION PATTERNS 
BY A VECTOR THEORY USING DIGITAL COMPUTERS 


INTRODUCTION 

Tile field equations 


E(x',y , z') 


1 1 

jojeo 


Vi/, H, • d 1 - [j (n X H, ) <// + (n • E, )V0 


dS 


and 




Jsi 


I ) xVv// dS 


are applicable to various geometrical and physical situations. In specific cases 
the Ej and Hj fields must be determined, usually from some radiation source or 
prime feed removed from the reflector surface. Expressions for the surface 
normal n and the differential area d S can ordinarily be found by standard 
mathematical techniques, but moderate amounts analysis are sometimes re- 
quired. The specification of the boundary edge d 1 is inherent in the definition 
of the surface. Only physical constants and other terms which are independent 
of the reflector geometry remain and integration can, in principle, be carried 
out by means of high-speed digital computers. 

Certain general observations can be made prior to the introduction of a 
specific reflector geometry, and the approach to the overall problem can be 
organized somewhat to facilitate computation Electromagnetic boundary condi- 
tions provide the transition from the and Hj fields to the incident fields E. 
and H. , An inspection of the field equations shows that the electric field appears 
in only one of the four integrals. If the reflecting surface is in the spherical 
wave zone of prime feed, which is the usual case, then it is possible to relate 
Ej directly to H. by a simple expression and cast the entire formulation in 
terms of the incident magnetic field H. . A consequence of this organization is 
the elimination of the unit electric polarization vector in the computations. It 
will be shown that the unit magnetic polarization vector is also simpler to ex- 
press and utilize than the corresponding electric vector. The selection of the 
magnetic vector is particularly significant since the field integrals will be treated 
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as summations of a large number* of weighted differential areas for every 
pattern field point. 


APPLICATION OF THE FIELD EQUATIONS 

The transitton from the fields E, and of the nonconducting medium to the 
incident fields E. and H. is made by regarding the antenna surface as locally 
plane so that the usual boundary conditions apply, 

n* ” 0 nx(E 2 -Ej) - 0 

CJj -0 

— ~ (7 ^--00 

n-(D 2 ~Dj) = T) 

The boundary conditions state that the normal components of B and the tangential 
components of E are continuous across the specified boundary. The normal 
components of D and the tangential comiDonents of H are discontinuous, however, 
for that boundary. Since all fields vanish within the perfect conductor, it follows 
that the incident £md reflected fields combine to satisfy the boundary conditions 
when B2 = H2 = D2 ^ E2 = 0 . Then 

* 

Hj-dl = 2H.-dT, 


nxHj = 2nxH. 
and 


n*Ej = 2n*E. 

I . . ■ • ‘ 

¥ 

in the field equations. * 

Electric and magnetic fields from the prime feed are related by the expres- 
sions^ 

IThe magnetic vector is calculated approximately 10^ times in computing the field intensity at one 
point of the radiation pattern of a large S.*band antenna. 

^Reference 1, page 284 
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and 



where is the direction of propagation at thej)rime feed. Here == |/mq /cq 

Ie| 

is the intrinsic impedance of free space, and ~ Zp. It is noted that the 

spherical -zone fields in the equations for Ei[x',y.', z') and H(x', yV,z') are orthogo- 
nal in space, synchronous in time, and also yield the intrinsic impedance of free 
space. 

The field equations can now be rewritten as 


E(x',y', z') r- - ^ ~ • dl - ^ J j (n xH. ).// 


( 1 ) 


~ {' t ) " • (a ’‘H j)Vv!']dS 


and 


H(x', y', z') 


f (nxH, 

Js. 


. )xVi/<dS, 


( 2 ) 


taking /ll^ = = 1 in the space in front of the reflector. At this point in the 

development the unit magnetic polarization vector is needed, and the history of 
the wave from the prime feed can be introduced. In addition, the history of the 
waves radiated from the reflector should be considered before proceeding to a 
specific reflector geometry. 

Assuming that the reflector is in the spherical wave zone of the source, the 
unit electric and magnetic polarization vectors e. and h. must be orthogonal to 
Pi and to each other. The unit electric vector will be assumed ^ to be coplanar 

^Ref. 2, page 3; Ref. 3, page 3 
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with an electric source polar isration vector v as shown (Figure 1). Then 

7T^x(vx“, ) 

since e. j.“, and c. i (vxT^). 


X 



For cases where magnetic fields are encountered, as in (nxH. ) of the field equa 
tions, it is necessary to write ^ 




at a given dS on Xj. 


*^Ref. 4, page 150 
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Effectively, 


‘ |pjx[p,x(vx^, )] I 

via ?. , 

1 

But 

Pi XV 

h. ^ r:~ir^ 

' ipjxv: 


directly by inspection of the physical picture or the reduction 


Ax [A X (B X A)] - “ A X [A X (A x B)] - A x B 

when A is a unit vector and B is any vector. The elimination of two of the 
three vector cross-products is achieved through the use of h. instead of e . for 
^e three field integrals which contain H. . The remaining field Integral contains 
E. (or pjxH. ) and two vector cross-products are necessary, using either e. or 
h. j to account for the polarization from the prime feed. The orientation of the 
electric source polarization vector v is arbitrary and it will be treated as a free 
vector^, not necessarily bound to origin 0, in subsequent developments. 

The history of the wave from the prime feed is now written into the problem 
explicitly. It is convenient to begin with the expression^ 


E. 


^Ref. 5 , page 35. in a sense v is **bound*’ to the feed at point /^in later developments. 
^Re/; 4, page 15C 
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where p is the distance fi’om feed to reflector, P is the total radiated power 
and Gf is the gain function of the feed. Capital letters 0, <I) are used for the 
spherical coordinate angles associated with the feed, and 0 , ip are reserved for 
the field points on a surface (usually a sphere) of observation. Then 


H. -* “ 


The l/p space divergence (attenuation) from feed to reflector surface and the 
accompanying phase change are now a part of the history of the waves 

designated E. and H. in the field equations. 


i^y 


4 ^ 




h. (0,4>) 


The history of waves from the surface of an antenna is contained in the 
scalar 0 and the vector of the field equations . Here . 


and 

Vl/; - (jk + -^) P 1^ . 

and r is the distance from the reflector surface to the point of the observer. 
These factors are under the integral sign for the general case. Distance r is 
variable for the variousi differential areas of surface , and vector 1^ is a free 
vector whose Cartesian components vary across as showia in Figure 2a. 
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If the observer is sufficiently far from the reflector , a spherical wave 
having only transverse field components results. This field decays as l/r and 
the phase delay from the surface > j to the field point can be written as e^i**^* 
to obtain the diffraction pattern. (See Figure 2b.) In this case the l/r divergence 
term can be taken outside of the integral sign, and the transverse components of 

Er(x , y , z ) - ~ r j a'/j,Q (n xHj ) 0d S 

(Equation 1, page 3) 

are obtained by forming . 

E^’ - E(x\y’,z') • 

and 




E(x',y',z') 




FEED TRANSLATION AND ROTATION 

The combined translation and rotation of a prime feed provides six additional 
degrees of freedom for controlling the illumination at a reflector surface. It is 
assumed that this surface has a parametric representation so that a translation 
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vector Cp^ ) and an orthogonal rotation matrix (A) can be introduced Jo obtaiii the 
correct magnitude, phase, and polarization of the comple'v vectors E. and H. at 
each dS of , Initially the source pattern function is taken to be F = 

(0,^>). See Figure 3. 
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and 


/:j X V p\ X V 

■' I —4 ■ M ■■■ -..i. 

X vl Ipj X V. 


The components of vectors in the coordinate frames with origins 0 and O' are 
indistinguishable and the indicated vector cross-product is well-defined. The 
intensity directed toward a point on surface can be found in terms of the dis- 
placement since F (©,<!>) - F (0',<^'). See Figure 3. Now 1 p' | = (x’^+ y'^ +z*2)^^ 
can be calculated, and the two unknowns ©',/*'' can be found from a system of three 
equations. 

X* = p' sin ©' cos 
y'. = p' sin ©' sin <!>' 
z' = p' Cos ©* 

If the tr^slation is degenerate, x^ - y^ = z^ == 0, and the angles ©' and <!>' 
become 0 and <^, respectively. Then F (0, 4>) is evaluated with the ordinary 
spherical coordinate angles, as before. 

If the feed is displaced and disoriented, or simply disoriented, the vector p^ 
and a 3 X 3 rotation matrix^ A with Euler angles a, y csui be used to determine 
the illumination at the reflector surface. Ordinarily the prime feed function F 
and the source polarization vector v are kno\vn in body coordinates, but the re- 
flector surface and the field points are designated in space coordinates. The 
inverse matrix transformation. 


r ^ 


■1 


r 1 

v" 

* 


(cos-ycosa-cos/Ssina siny) (-sinycosa-cosPsimcosy) (sin/? sina) 


' v; ■ 

'^y 

- 

A-' 


v" 

y 

= 

(cosy sina+eos/?cosa sin>) (~sinry sina+cos/?cosacosy) (~sin/3cosa) 


v" 

y 

\ 






(sin/?siny) (sin/3cosy) (cos/?) 


v" 

z 


space body 


^Ref. 6, page 107 
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provides the components of the source polarization vector in space coordinates 
so that the magnetic polarization vector becomes 






p, XV 


\p. X V 


1 


The vector cross-product is well-defined when the components of p[ and v“ are 
all relative to space coordinates. 

The intensity of the feed pattern F which is directed toward the reflector at 
a point (x,y,z) in space coordinates can be determined by the direct rather than 
the inverse transformation. Vector p' has components which are identical in 
both 0 and O' coordinate systems. The components of p' relative to lx)dy coor- 
dinates, frame 0" , can be found via the matrix transformation 




■ 



v 


y” 

- 

A 


y’ 

¥ j 



— 


_z'^ 

where 


A-i 

— • 

AT 



since the rotation matrix is an orthogonal matrix. 

But 

x" = sin0” cos4)" 
y" - p" sin0" sin$" 
z" ^ p" COS0” 

again provides a system of three equations in two unknowns, 0", and If the 
rotation is degenerate ( a = /3 =7 = 0) so that A = I , the identity matrix, the 
equations for the case of simple feed translation reappear. 
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The preceding methods for introducing feed translation and rotation into the 
problem are adapted to numerical computation rather than analysis. The re- 
flector area is presumed to be subdivicTed suitably into a number of differential 
areas. If the subdivision or sampling^ is adequate, then the field intensity due 
to the illumination can be used to ’’weight” each differential area prior to sum- 
mation over the reflector. The precise manner^ in which the area is ’’weighted” 
is a topic sufficientlj^ broad to deserve detailed treatment, and will not be con- 
sidered in this document. It remains to provide the general methods for obtaining 
the local normals (n) and the differential areas (dS) on reflector surfaces. The 
rotationally symmetric surfaces— paraboloid, hyperboloid, ellipsoid, sphere and 
cone— are presented in Appendix A as they are representative cases which find 
frequent application in reflector antenna systems. 


SUMMARY 

This report provides mathematical subroutines which are required tb trans- 
port the illumination of a prime feed to a reflector surface. The feed translation 
vector and the rotation matrix afford an unusual amount of flexibility to the pro- 
gram. Utilization of the magnetic quantities throughout greatly simplifies the 
formulation and external control. Although the appendix contains only subroutines 
for the normal and differential area of rotationally symmetric surfaces derived 
from the conic sections, the field equations are not restricted to these particular 
surfaces. Other smooth surfaces, symmetric. or non-symmetric, can be intro- 
duced providing that their sizes and minimum radii of curvature exceed one- 
wavelength. 

Dual or multiple-reflector systems can be analyzed by the techniques outlined 
herein. Complex prime-feed arrangements can also be accommodated by super- 
imposing field solutions or. prime feed-functions. In conclusion, a library of feed 
arid reflector-surface subroutines can be formed and utilized for many practical 
configurations. 
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^Ref. 7, page 323 
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APPENDIX A 

Normals and Differ ential Areas from Parametric Equations 
^Cylindrical Variables (cTj t^) 

X = a- sin I • y = - c 7 cos ^ 

paraboloids , . d , 


hyperboloids 

ellipsoids 

spheres 

cones 


Z c(l + cr^/a^)’^ + Zj 
z = c(l - +■ Zj 

2 = (c^ “ + Z^ 

Z = ccr + Zj 


iMgl^nts to Surface : = ||. p(x,y,z) = i x + Jy + kz = (x.y.z) 

= (cr cos cr sin 0) 


paraboloids 

= (sin - cos cr/2F) 

hyperboloids 

p^ = [sin - cos , crc/a(a^ + 

ellipsoids 

p^ = [sin -'COS - crc/a(a^ - 

spheres 

py = [sin - cos C. ~ cr/(c2 - 0 - 2 )'^] 

cones 

p^ = (sin ^ , ~ cos 5, c) 

Normal to Surface: 

rT = (p^ X p^)/|p^ X p^ 1 

Differential Area: 


dS = (EG - F)« do- d( 


E = X / + .yj + il, F ^ 

= x^Xr+y^y^ + z^z^, G = + y^2 + ^2 

paraboloids 

dS = C7'(l + (j^/ 4F^)'^ dcr d{, 
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hyperboloids 

dS = 

!7 [l + 

ellipsoids 

dS = 

O' [l + cr2 c2/a2(a^ ~ cr^)] 

spheres 

dS = 

c7c/(c^ - CT^)’^ da dl 
1 

cones 

dS 

a(c^ + 1)*'^, da d(. 

See also: 




Ref. 8, page 3 
Ref. 5, page 206 
Ref. 9, page 3.5-10 
Ref. 10, page 106 







” dcT dT 
^ dcr 
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APPENDIX B 

Minimum Radius of Curvature from Parametric Equations 
in Cylindrical Variables (cr, { ) 


Curvature k, = X/R , where R = radius of curvature.' 

l C ' c 

X ^ X (cr) , y = y (cr) , z - 7 . (cr) 



2 

.H _ 

2 

•< 

^ y, 


y z 


Z X 


+ 


+ 


1 1 


1 1 


1 ■<! 

X y 


y z 


Z X 


(1 + Z' 2)3 


(X 


• 2 


.• 2 


+ z’ 2)3 


Since 


- sin ^ ♦ x" 0 , y‘ - ~ cos ^ , y" = 0 , z' = k • 


kj - maximum at o- - 0. 


paraboloids 


z" = 1/2F, 


min R. 


hyperboloids 


•z - 


a 




[(a2+cr2)3/2 (a?+c72)i/2j 


min R. 


ellipsoids 


z - ~ 


+cr- 


1 


(a2-c72)3/2 ^g2_^2^1/2 


— min R^ 
2 c 


spheres 


z” = 


“0-2 ' 


1 


(c2-^2)3/2 (^2«^2)l/2 


min R. 


cones 


z" = 0 


min R^ 


'Ref. 10^ page 51 
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= 2F 



c 


- c 


= 00 



